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I ^i ' Abstract. We compute the probability of finding metastable states at a given field 

in the mean- field random field Ising model at T = 0. Remarkably, this probability 
is finite in the thermodynamic limit, even on the so-called "unstable" branch of the 
magnetization curve. This implies that the branch is reachable when the magnetization 
is controlled instead of the magnetic field, in contrast with the situation in the pure 

^D ' system. 

cn- 

o 

OO 

O 

> 



Stable, metastable and unstable states in the mean-field RFIM atT = 2 

1. Introduction 

As is well known, the main approximation underlying mean-field theories of phase 
transitions consists in neglecting any spatial dependence of the order parameter. In 
simple systems such as fluids or ferromagnets this leads to an equation of state 
(e.g. the van der Waals equation) that exhibits a continuous loop below the critical 
temperature T^. In this framework, one can distinguish between stable (equilibrium), 
metastable, and unstable states. The intermediate, unstable branch of the loop is 
associated with maxima of the free energy and has a negative slope (whence a negative 
susceptibility), and the unstable and metastable regions of the phase diagram are 
separated by a spinodal line where the free-energy barrier vanishes. Whereas the concept 
of metastability can be in some sense extended to systems with short-range interaction 
in finite dimensions where it becomes a matter of timescales, the unstable branch is 
a complete artifact of the mean-field approximation and the homogeneous states along 
this branch do not represent any real physical situation. 

Things are different in the presence of disorder since inhomogeneities are induced 
even at a microscopic scale. As a result, configurations formed by a multitude of small 
domains of the two phases can become metastable in an extended range of the "external" 
field (e.g. the fluid pressure or the magnetic field). Interestingly, this can already occur 
at the mean-field level and the purpose of this paper is to emphazise that metastable 
states are present all along the so-called unstable branch in the mean-field random 
field Ising model (RFIM) at T = 0. This feature is in fact a precursor of what happens 
in the RFIM in finite dimensions, and it has consequences for experiments in actual 
random- field systems at low temperature. 

2. The T = RFIM on a fully-connected lattice 



We consider a collection of A^ Ising spins (sj = ±1) interacting via the Hamiltonian 

J 

'2N 



^ = -^E^^^^-E(^^+^)^^ (1) 



i^j i 

where J > 0, if is a uniform external field, and {h^} is a collection of random fields 
drawn identically and independently from some probability distribution V{h). (In the 
following, we consider a Gaussian distribution with zero mean and standard deviation 
A.) This mean-field model can be obtained from the usual short-range RFIM by placing 
the spins on a fully-connected lattice (the connectivity is then equal to A^ — 1 and the 
exchange interaction is rescaled by A^ to ensure a proper thermodynamic limit). Eq. 
([1]) can be also rewritten as 

H = - 5^ f^S^ (2) 

i 

where fi = J{m — Si/N) + hi + H is the effective field acting on each spin i and 
m = {Y2- Si)/N is the magnetization per spin. 
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It is straightforward to compute the equihbrium properties of this model by the 
rephca method, without all the complications that plague the case of random exchange. 
In particular, the average magnetization in the thermodynamic limit A^ — >■ oo is solution 
of the self-consistent equation [1] 



m{H)= f V{h)tanh{p{Jm{H) + H + h)]dh 



(3) 



where (3 = 1/kBT. (In Ref.[l], only the case if = is considered, but the generalization 
to if 7^ is straightforward.) At T = 0, the equation becomes 



mo(ii) 



fvih) sgn{Jmo{H) + H + h)dh 
2p{mo{H)) - 1 



(4) 



where sgn(x) = x/\x\ and p{m) = jTu_j^V{h)dh. This yields mo{H) = erf([ii + 
Jmo{H)] /A^/2) in the case of the Gaussian distribution, where erf(a;) is the error 
function. Since the magnetization per spin is a self-averaging quantity, mo (if) is also 
the expected value of the magnetization in a very large sample. 
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Figure 1. Metastable states in tlie mean- field Gaussian RFIM for a single disorder 
realization of size iV = 5000 with A = 0.5 (the increment in the field is AH = 25.10""'). 
The dashed curve represents the solution of Eq. (j4|) . The inset shows that there may 
exist a few states with different magnetizations at the same field. 



Below the critical disorder A° (A° = ^/2/ttJ for the Gaussian distribution), 
Eq. (jlj) has three solutions in a certain range of the field and the curve mo (if) 
exhibits a characteristic "van der Waals" loop with an intermediate branch along which 
dm/dH < (see Fig. 1; we take J = 1 in all figures). As usual, this behavior is 
associated with the non-convexity of the (free) energy. At a given field H, the ground- 
state corresponds to the solution with the lowest overall energy whereas the intermediate 



Stable, metastable and unstable states in the mean-field RFIM atT = 4 

branch has to the largest energy. In the space of rephca magnetizations [T], such branch 
corresponds to an absolute maximum. The result from the replica method can be shown 
to be exact [1], and it seems natural to describe the intermediate branch as "unstable", 
naively transposing the situation found in the pure system (i.e. in the absence of random 
field). However, this is misleading. Indeed, since Eq. (jlj) results from an average over 
disorder (which restores translational invariance) when A^ — > oo, it says nothing about 
the orientation of the spins in a given (finite N) sample. It turns out that some of the 
states along the intermediate branch are metastable (i.e. local minima of the energy) 
and not unstable. Indeed, Eq. (j4]) is trivially verified if 

s, = sgn(/,),2 = l...iV, (5) 

which is the definition of the so-called one-spin-fiip stable states. There is nothing new 
here: it is known that Eq. (jl]) also describes the nonequilibrium (hysteretic) behavior 
of the mean-field RFIM at T = with the single-spin-flip (Glauber) dynamics[2]. As 
H is slowly varied from ±oo, this dynamics imposes at any field that the spins with 
hi < —Jm — H point down whereas the other ones point up, and the self-consistent 
equation for the average magnetization, m{H) = JV{h)sidh, is just Eq.(|l]). (Note 
incidentally that this corresponds to an "annealed" average, but it yields the correct 
result in this case.) The nonequilibrium RFIM at T = has been the subject of extensive 
studies in recent years [3] and the mean- field model, despite some peculiar features (the 
equilibrium and nonequilibrium critical disorders coincide and there is no hysteresis for 
A > A°), has the advantage that many interesting properties can be computed exactly 
(for instance the size and duration of avalanches) [21 H] . For A < A°, the nonequilibrium 
system explores the lower (or upper) branch of the loop until H reaches one of the two 
coercive fields where dm/dH diverges and the magnetization jumps discontinuously. 
This "infinite avalanche" [2| occurs when each flipping spin triggers one other spin on 
average, which corresponds to the condition 2JV{—JmQ{H) — H) = 1. Hereafter, we 
shall use the shorthand Vq for the quantity V{—Jmo{H) — H). 

The point that we want to emphasize here is that a configuration where each spin 
satisfies Eq. ([5]) is a metastable state by definition, whatever the value taken by mQ{H) 
(of course, a state on the intermediate branch has a larger energy than the corresponding 
ground state and it cannot be reached by controlling the field). This is illustrated in 
Fig. 1 that shows the magnetizations of the metastable states in a single sample of 
size N = 5000 with A = 0.5 as a function of H. One can see that there are stable 
states in the intermediate region and that they gather along a curve that will become 
the so-called "unstable" branch in the thermodynamic limit. At a given field H, the 
number of these states is very small, obviously not exponentially growing with system 
size. However, the probability of finding metastable states remains finite when N — >■ cxd, 
as we now show. 
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3. Metastable states 

For a disorder realization of size N, consider the ensemble of metastable configurations 
at the field H with exactly P spins up (and thus an overall magnetization M = 2P — N). 
The number of such configurations is 

Af{M, H) = Tr^.^i W Q{siM)6K{Yl ^i - M) , (6) 

i i 

where 6(x) is the Heaviside step function and 5k is the Kronecker-5. It is easy to see 
that there cannot be more than one metastable state with magnetization M in a single 
realization (i.e. M{M,H) = or 1). Averaging over disorder, we have 

U{M,H) = Cp)p{m - 1/NY[1 - Pirn + l/iV)]^"^ , (7) 

where the two terms in the right-hand side represent the probabilities of having P spins 
up and N — P spins down at the field H, respectively. Using the Stirling approximation 
for the factorial and expanding p{m ± 1/A^) to first order in 1/A^, we find 



with 



JfiMjT) ~ . /X4^!!^e-^^(-^-^'")[^+^TS^] (8) 

V ttN y/l - m? 

,, , 1 + m, 2p{m) 1-m, 2\l - pirn)] .. 

(him) = n ^ ^ H In — Ll^^l . (g) 

^^^ 2 1 + m 2 l-ni ^ ' 

This number is exponentially small {(pim) < 0) except when m = 2p{m) — 1 = molH), 

i.e. when m is solution of Eq.(jl]): (f){m) is then maximum and equal to 0. Expanding 

(/)(m) to second order close to m = mo{H), we find that the average number of metastable 

states at the field H is finite when N ^ oo and is given by 



^f(H) = J2^iM, H)r^— dm M{M, H) 

M -^ 



|l-2JPo1 



(10) 



This result is valid both above and below A°, and for A < A° the three branches of 
mo{H) must be considered separately (with 2JVq > 1 along the intermediate branch). 
The behavior of M{H) as a function of mo{H) is shown in Fig. 2 for A/ J = 0.5. It is 
worth noting that Mi^H) diverges at the coercive fields (i.e. at the spinodal endpoints) 
where 2p'{m) = 2JVq = 1 and 0"(mo) = -[1 - 2JPo*]V(l - mg) = (from now on, 
the dependence of tuq on H will not be indicated for brevity). In this case, the large- 
deviation function (f){m) must be expanded to fourth order about m = itlq and N'{H) 
scales like N^/^. At the critical point (A = A° and H = mo = 0), 0(m) must be 
expanded to the sixth order and N'iH) scales like N^^^. The special form of 0(?7i) is 
responsible for these unusual mean-field exponents. 

Note that it is crucial to take into account the contributions of order 1/A^ in the 
local fields /, in order to obtain the correct result. This is due to the fact that the 
function 0(m) is zero at the saddle-point, so that the associated complexity is zero. 
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Figure 2. Average number of 1-spin-flip (solid line) and 2-spin-flip (daslied line) stable 
states along the curve mQ{H) for A = 0.5. Both quantities diverge at the spinodal 
endpoints. 



(This is different from tlie situation found in tlie mean-field spin glass model discussed 
in Ref. [5].) On tfie otfier fiand, tfie subdominant terms play no role in determining the 
probability of finding a metastable state with magnetization m at the field H, which is 
only given as usual by the fluctuations around the saddle point: 



Af(m, H) 



N(j)"{mo) ^N^"(mo)(m-mo)^ 
6 ^ 

27r 



:iii 



oo if 0"(mo) ^ 0. 



AfiH) 
This becomes a (5-distribution when A^ 

More generally, we can compute all the moments M{H)^ and use this information 
to obtain the full probability distribution P{q) = 5k{M{H) — q) (with J\f{H)^ = g" = 
^g"P(g)). This amounts to count all possible ways of ordering n magnetizations 
Mi,M2...M„, given that J\f{Mi, H)M{Mj, H) = J\f{Mi,H) when M, 
7V(M, H) is just or 1. This yields 



Mj since 



A^(^)" = E E in,,n,...nr)\J2l 



=1 ni,n2,...nr>l 



Afi,M2...,M, 



Af{M,, H)Af{M2, H)...X{Mr, H) (12) 



where (ni, n2---nr)\ = n\/{ni\n2\---nr\) is a multinomial coefficient and the sum runs over 
all (ni,n2, ...nj.) such that Yll=i''^i ~ ^ ^^^ ni,n2...nr > 1- The notation ^^ indicates 
that the sum over the magnetizations Mi, M2..., M^ is restricted to a specific order, say 
Mr > Mr-i > ...M2 > Ml. We then define Mi = 2Pi - N, where Pi is the number of 
spins up, and introduce the (strictly) positive quantities Qi = Pi — Pj_i (z = 2, ...r). It 
is easy to see that Eq. (^^ generalizes to 



Af{Mi, H)Af{M2, H)...Af{Mr, H) = N\- 



Pil 
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r 



\p{mi - 1/N) - p{mi_i + l/iV)]Q» [1 - p{mr + 1/N)] 



n-Pi-J:U2Q 



fj QJ (iv-Pi-EL.Q.)! • ^ ^ 

Since only the magnetizations in the close neighborhood of ttiq contributes to the sum 
E^ when A^ ^ 1 (with the Qj's being at most of the order N°' with < a < 1), we can 
expand p(mi±l/A^) aroundp(mo) to first order in 1/A^ and use 1 — Pi/N ~ (1— mo)/2 = 
1 —p{mo). After some straightforward manipulations, we then obtain 

iToo Em,,m....,m. mMuH)AriM,,H)...Ar{M^,H) = ATiH) 

x\j2^^^i^J'Poe-''^^'fr' (14) 

Q>i ^' 

and, after inserting this expression into Eq. (1121) . 



^^^ = E E ini,n2,...nr)\[a{mo)-lY-' (15) 

-'^(-nj r=l ni,n2,...nr>l 

where 

«(^o) = E ^^^(2JPoe-^^^")'^ . (16) 

fc>0 

We recognize in ( iT6l) the series expansion of 2;/[iy(2;)(l + iy(2;))] near the origin, 
where W{z) is the so-called Lambert function, defined as the root of the equation 
W{z)e^^^^ = z^. This series converges for \z\ < 1/e, which is always true in Eq. flT6l) 
where z = — 277^0 e'^'^^o. (The series f lTB]) only refers to the principal branch 1^0(2^) 
which takes on values between —1 to +00 for z > —1/e and is analytic at z = 0.) As a 
result, 

-2JP*e-^-^^o 



which yields 



e-2j-Po* 



a{mo) = Y^l^JpZ = ^{H) if 2 JPo* < 1 ^ 



«("^o) = , , T,, , o7^*.-2jPn if 2 JP* > 1 



1 + Wo{-2JV^e 

Knowing from Eq. (IT5l) all the moments g" of the probability distribution P(g), we 
can build the generating function 

5^e-M?) = i + E(-^r^" (19) 

g>0 n>l 

to obtain 
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Figure 3. Probability of finding q metastable states in zero field for A = 1 (blue 
circles) and A — 0.5 (red squares) in the mean-field RFIM (the statistics is taken over 
5000 disorder realizations of size N = 20000). For A = 0.5, only the states in the 
vicinity of the intermediate branch (around m — 0) are counted. The dashed lines are 
guides for the eye. (Color on line) 



with q = N'{H). This equation can be inverted, showing that P{q) decreases 
exponentially for g > 1. More precisely, we have 



P(0) 
P{q) 



a 
Q 



for g > 1 



(21) 



a{a — 1) a 
where both q and a are functions of mQ{H). 

For A > A°, a = g so that P{0) = and the most probable value of Mi^H) is 
g = 1, as could be expected. (What is perhaps less expected[7] is that P{q) 7^ for 
g > 1 and that g = M{H) > 1.) For A < A|?, the most probable value of M{H) on 
the intermediate branch is g = except very close to the spinodal end-points where 
it is again g = 1 (for 1 < 2JP* < 1.0073, P(l) = q/a^ > P(0)). Note also that 
P{q) decreases more and more slowly when approaching the spinodals as the inverse 
characteristic scale ^~^ = ln(l — 1/a) — > 0. In all cases, there is a finite probability of 
finding a few metastable states at a given field^, as illustrated in Fig. 3 that results 
from an exact enumeration of all metastable states in 5000 disorder realizations of size 
A^ = 20000 at H = for A = 1 and A = 0.5 (in the latter case, only the states in 
the vicinity of the intermediate branch, i.e. around niQ = 0, are counted). One has 
a = q ^ 2.12 for A = 1, and g ^ 0.340, a ^ 1.324 for A = 0.5. The numerical data 
shown in the Figure are in very good agreement with the predictions of Eq. fl2T]) . 

The above calculations can be generalized to 2,3, ...fc-spin-fiip stable states, i.e. to 
spin configurations whose energy cannot be lowered by the flip of any subset of 1, 2, .../c 
spins ^. It is easy to see that a conflguration with P spins up and N — P spins down 
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is A;-stable if it is {k — l)-stable and if the random fields on the P spins up satisfy 
X]a=i ^ic > —k{Jm + H) + k'^J/N whereas the fields on the N — P spins down satisfy 
Sq=i hi^ < —k{Jni + H) — k'^J/N for any subset {ii,i2, •••4}- From this, one can for 
instance compute the average number of 2-spin-flip stable states at the field H and find 
that 

-2JV, 



[2e-^jvs _ g- 



^^^^(^)- ^ |1-27P,1 ^ ^''^ 



when A^ ^ oo. The comparison with Af'^'^^H) = M{H) is shown in Fig. 2. Of course, 
one has N'^'^\H) < N''^'^^^\H)... < J\f^^\H). We have not investigated the behavior for 
k ~ viV|10j. In any case, in order to go from a metastable state on the intermediate 
branch to the ground state on the lower or upper branch, one needs to flip a number 
of spins of order A^ (which corresponds to take into account another solution of the 
saddle-point equation). 

4. Conclusion 

In this paper, we have computed the probability of finding metastable states along 
the so-called "unstable" branch of the mean-field RFIM at T = and shown that the 
average number remains finite in the thermodynamic limit. 

The presence of a few metastable states along the intermediate part of the curve 
mo(i?) for A < A^ can be considered as a precursor of the phenomenology observed in 
the T = RFIM with finite-range exchange interaction. This will be investigated in 
a forthcoming paper fTT] dealing with random graphs of large but finite connectivity 
z. Preliminary results indicate that, as soon as z is finite, a strip of finite width 
develops around the curve mo{H) in the field-magnetization plane, strip in which the 
density of the typical metastable states scales exponentially with the system size. This 
occurs both above and below A^. As the connectivity decreases, one expects the strip 
to widen but to remain distinct from the actual hysteresis loop in the weak-disorder 
regime. In fact, as suggested in Refs. [T2| [13] . one may associate the discontinuity in the 
hysteresis loop below A^ to the existence of a gap in the magnetization of the metastable 
states beyond a certain value of the field. Of course, all this is strictly valid only at 
T = 0. However, as is well known, free-energy barriers are very large in random- field 
systems and thermally activated processes are not expexted to play a significant role 
on experimental time scales, at least at low temperature. Therefore, the above picture 
is expected to be relevant to real situations. In particular, the presence of metastable 
states (and not simply unstable ones as in pure systems) in the central part of the 
hysteresis loop means that this region could be experimentally accessible, for instance 
by controlling the magnetization instead of the magnetic field (and more generally the 
extensive variable conjugated to the external field). This can be put in relation with 
the re-entrant hysteresis loops that are observed in some magnetic systems[Tl] or in 
shape-memory alloys [15] (see also the discussion in Ref . [I6] ) . 
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